| 


“ AD=A063 537 ARMY ARMAMENT RESEARCH AND DEVELOPMENT COMMAND ABERD=-ETC F/6 18/3 
USER'S MANUAL FOR CASSANDRA: CLOUD SNAPSHOTS OF DUST RAISED ALO--ETC(U) 


NOV 78 RL SHOWERSs C CRISCO 
UNCLASSIFIED ARBRL-TR-02116 SBIE-AD-E430 152 NL 


FILE COPY 


a 


t ee: 
,O2 
ics 


UNCLASSIFIED 


SECURITY CLASSIFICATION OF THIS PAGE (When Date Entered) 


RE. s ONS 
REPORT DOCUMENTATION PAGE 
. REPORT NUMBER / Ee iy GOVT ACCESSION NO. 3. RECIPIENT’S CATALOG NUMBER 


TECHNICAL REPOR 
5. TYPE OF REPORT & PERIOD COVERED 


4. TITLE (and Subtitle) 


USER'S MANUAL FOR CASSANDRA: CLOUD 
SNAPSHOTS OF DUST RAISED ALOFT 


| Final | ‘ 


6. PERFORMING ORG. REPORT NUMBER 


8. CONTRACT OR GRANT NUMBER(a) 


H 
i 
' 


ite AUTHOR(a) 
\casi’cs L. /Showers |; 


jAB, APES / 


i 


Carl ‘Crisco } 


10. PROGRAM ELEMENT, PROJECT, TASK 
AREA & WORK UNIT NUMBERS 


9. PERFORMING ORGANIZATION NAME AND ADDRESS 


U. S. Army Ballistic Research Laboratory 


ATTN: DRDAR-BLV RDTGE|1W162118AH75 Vv “ALY 
Aberdeen Proving Ground, MD 21005 Lt a : 


11. CONTROLLING OFFICE NAME AND ADDRESS AZ. REPORT DATE _ 
U. S. Army Armament Research § Development Commarid|, 
U. S. Army Ballistic Research Laboratory ( 
pile DRDAR- BL 
p D 0 


15. SECURITY CLASS-e ia feport) . ” 


| UNCLASSIFIED 


'Se. DECLASSIFICATION/ DOWNGRADING 
SCHEOULE = 
Ke ~ of 


Approvea for public release; distribution unlimited. ih oafe nn 


16. DISTRIBUTION STATEMENT (of thie Report) 


17. DISTRIBUTION STATEMENT (of the abstract entered in Block 20, if different from Report) 


18. SUPPLEMENTARY NOTES 
This research was jointly sponsored by: (a) AMC-Army Project 1W162118AH75, 
Task AB, Nuclear Weapons Environments, Phenomenology, and Effects Studies- 


Fallout, (b) Defense Nuclear Agency-under Subtask No. V99QAXNAO11, Work Unit 
No. 01, Residual Radiation Environment - Fallout Definition. 
19. KEY WORDS (Continue on reveree side if necessary and identify by block number) 


Radioactive Fallout Fallout Exposure 
Collateral Damage Dust Clouds 
Fallout Prediction 


20. ABSTRACT (Continue em reverse side If neceeeary and identify by block number) 
‘The computer code, CASSANDRA)|has been developed to simulate the spatial distri- 
bution and concentrations of dust raised aloft by a nuclear cloud. The code 
utilizes the DELFIC cloud rise description as the basis for a dynamic wafer 
shape/transport model. Dust concentration for a particular particle size class 
is calculated directly from the radius of the contributing wafer at the altitude 
of interest. 

The outputs of the program are either an array of numbers specifying the values 
of collective dust densit gridwork of points in space or a tabulation of 


DD anes 1473 * eornow oF 1 wov 6515 OBSOLETE UNCLASSIFIED 


— 
SECURITY CLASSIFICATION OF THIS PAGE (When Data Entered) 


UNCLASSIFIED 
én, -/ / SECURITY CLASSIFICATION OF THIS PAGE(When Date Entered) 


dust concentration at specified coordinates and times. 

As a test case a simulation of the high explosive test Dial Pack was perfotned. 
The computed values were within the range of values reported from the cloud 
sampling. 

In addition to a general description of the code, a user's guide has been 
incorporated. Examples have been included of computer generated output along 
with a complete listing of the CASSANDRA code. 


| 


\ 
i\ 
N 


UNCLASSIFIED 


Ce ee ee 
SECURITY CLASSIFICATION OF THIS PAGE(When Date Entered) 


TABLE OF CONTENTS 


Page 

LIST OF ILLUSTRATIONS .. .... ccc cccccccccerccccreccccncccees ah ed) 

I, INTRODUCTION.......---.00. sahatole cioieiaistsialelete cielsialele aver eieioiale/clatecsteleie 
II. CODE DESCRIPTION. ... cc ccceeceeecvcvccceces A OBS BOGOR OS 5 
III. FORTRAN LISTINGS OF CASSANDRA......... aiel/atalersiiai#i0'sl'arefeteie ajatscotsier (2 
IV. USER'S GUIDE TO CASSANDRA........+-6- Sieieleerere eels sie aeislelats siete erste 0 
A. Control Cards for UNIVAC 1108 EXEC 8........eeeeee- cries 130 

B. Data Deck for Execution of CASSANDRA.........eeeeees waters el OU 

V. EVALUATION OF CASSANDRA........ rel aiehsielas senna ale sialereieleleiereene ceks « «L358 
A. Dial Pack Simulation, ......cccccccccccccccrcsccsccscccens 138 

B. Sample Data Deck for Dial Pack...........0-. are rareie ate sieie es 141 

C. Sample Computer Output for Dial Pack..... AO BAEOR SO EOS ONES 142 

DB. Diak Pack<Simulation Results. oc. .<ccss.c cose cos cess o0cv kl 
DISTREBULION (LEST iiiiscce ctecieisisisisie scien ss AGO COUOCE DEG COICO ORE 167 


LIST OF ILLUSTRATIONS 
Figure Page 


1 Cross Section of CASSANDRA Cloud 1.0 kt at 2.0 Seconds......... 7 


2 Dust Concentration at Target Coordinates at Specified 


SUDO terarat ate) ec oraiaycrarsie er siere cia cree neal ete Bieieiareleisiaisisis sjavialejlelelein/s\els elie isiejee\e e@ 9 
3, Dust Concentration at 350 Seconds < sci s.5cis1es « nejs om seis sin's + aj09 0 seu 
A CUOUS RUS Cr rvehevateheie erelerereeresie Bier eterete Sielieie elise )e se | ¢ 01 ele 0160/0, e/\ejs\siecns ee 139 
S Dial Pack Particle Size Distribution... s0cci0escaecces o's 0c el40 
6 Dial Pack at 615 Seconds Dust Concentration.....cesececscceeces 162 
7 Dial Pack at 804 Seconds Dust Concentration........ceeeereeee +0163 
8 Dial Pack at 950 Seconds Dust Concentration............- we supe 0 164 
9 Dial Pack at 1304 Seconds Dust Concentration. .......seeeee% 2022165 


10 Dial Pack at 1500 Seconds Dust Concentration............ os 0% se 3166 


I. INTRODUCTION 


A nuclear detonation on or close to land creates large persistent 
clouds of particulate debris, or dust, consisting of dense particles 
with sizes ranging to thousands of microns. Such a debris cloud poses 
a threat to RV's reentering the atmosphere at high velocities or ABM's 
leaving the ground. To determine the threat of erosion of missile heat 
shields, one must know the spatial distribution and concentrations of the 
Material in the air. The computer code CASSANDRA has been developed as 
a means of supplying that information. 


CASSANDRA is an acronym for "Cloud Snapshots of Dust Raised Aloft". 
It is a computer program designed to simulate the movement of dust 
raised aloft by a nuclear cloud. It offers a rapidly executed calcula- 
tion of dust distribution and concentration for any time after the cloud 
begins to rise. 


The code incorporates the highly detailed rising cloud model of 
the DELFIC! (Defense Land Fallout Interpretive Code) which provides a 
table of cloud parameters at predetermined time intervals from the time 
the cloud begins to rise to the time of stabilization. Other than the 
code necessary to produce the table of cloud parameters, CASSANDRA uses 
the particle fall procedures and atmosphere model of the DELFIC, but 
the remainder is a complete departure from that code. 


II. CODE DESCRIPTION 


The analytical and geometric basis of CASSANDRA begins with the 
description of the dust cloud at the initial time, which is the time 
that it begins to rise-usually several seconds after burst. For the 
purposes of cloud-rise calculations, the cloud is assumed to be an 
oblate spheroid. For dust calculations, the cloud is modeled, at the 
initial time, as a right circular cylinder. This cylinder is loaded 
with an uniformly distributed soil mass which is assumed to be composed 
of particles with an assigned size distribution. The soil mass in the 
cloud is subdivided into, at most, 200 particle size classes of equal 
mass and each of the particle size classes defines a separate, geometric- 
ally identical, colocated, dust cloud. Next, a geometric subdivision 


1. "Department of Defense Land Fallout Prediction System", DASA 1800 
a. Voi I - System Desertption, 27 Jun 66 (AD 483 897) 
b. Vol II - Initial Conditions, 30 Sep 66 (AD 803 144 
ce. Vol III - Cloud Rise. 19 May 67 (AD 819 770) 
Vol III - Cloud Rise, Revised, 1 Sep 70 (AD 879 890) 
d. Vol IV - Atmospherte Transport, 2 Feb 67 (AD 815 263) 
e. Vol V = Particle Activity, Feb 68 (AD 882 239) 
f. Vol VI - Output Processor, 20 Feb 67 (AD 814 055L) 
g- Vol VII - Operator's Manual, Avr 68 (AD 836 871) 


is made of each of the colocated clouds by passing through them a set of 
equally spaced, parallel, horizontal planes which divide the clouds into 
a series of smaller subcylinders which are referred to as wafers. The 
wafers are the basic units of transport of CASSANDRA. Each of the 
geometrical wafers, then, becomes a set of superimposed wafers, one for 
each particle size class, and after the initial time each of the wafers 
is assumed to be transported independently. For the purpose of defining 
a dynamic wafer shape, each of the wafers is structured at the initial 
time by an arbitrary number of equally spaced altitude reference planes 
parallel to the top and bottom. The altitude reference planes subdivide 
the wafer vertically into compartments which aid in defining an outer 
shape for the wafer as it falls from the cloud. 


Figure 1 represents, in cross section, the cloud at the initial time 
from a 1 kt detonation. The cloud has been subdivided into three wafers. 
The lower horizontal line represents the ground and the cloud has these 
dimensions and altitude at two seconds (initial time). The horizontal 
dashed lines represent an arbitrary number of internal altitude reference 
levels assigned to each wafer. These altitude reference levels provide 
internal reference altitudes at which wafer radii are calculated as 
different levels fall below the rising cloud bottom. These reference 
levels allow the wafer to assume a dynamically determined shape derived 
from the dynamics of the rising, expanding, wind-shifted cloud. 


The growth, shape, and transport of wafers can best be visualized 
by first considering the simpler case of a single dust particle in the 
cloud at initial time. As the gas cloud rises it will impart an upward 
velocity to the particle which is equal to the velocity of the cloud it- 
self. At the same time, the particle will have a downward component of 
velocity due to gravitational settling. The magnitude of the gravita- 
tional settling rate will be determined by the size of the particle, its 
mass and the density of the gas in the cloud at the time of interest. 
The particle, then, at any time, will have a vertical velocity which is 
the resultant of its upward and downward components. If the magnitude 
of the downward component is great enough, the particle will fall below 
the cloud bottom sometime before stabilization. 


The wafers are transported by assuming that the wafer tops and 
bottoms are independent and that fall rates are determined by assigning 
the largest particle in the particle size class, characterizing the 
wafer, to the bottom and the smallest particle to the top. The altitudes 
of the wafer top and bottom are calculated for the end of a particular 
time step by calculating rates of fall of the appropriately sized 
particles within the cloud and subtracting those from the rate of rise 
of the cloud. If the bottom of the wafer does not fall below the 
bottom of the cloud during the time step, then new altitudes of the 
intermediate reference planes are calculated by allowing them to be 
equally spaced between the new wafer top and bottom altitudes. Addition- 
ally, the wafer radius at the top and bottom and at each reference 
altitude and horizontal coordinates of the wafer top and bottom is 
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assigned the same as the cloud radius and horizontal coordinates at the 

end of the time step. If, during any time step, the wafer bottom falls 

below the cloud bottom, then the assumption that all variables are 

linear during the time step allows an exit time to be calculated. This 

time is used to calculate the cloud radius and horizontal coordinates 

which are assigned to the wafer bottom. The wafer bottom radius after 

exit is assumed to be fixed and the wind transport coding calculates the 

horizontal shift of the wafer bottom from the time of cloud exit to the 

end of the time step. Radii for the intermediate reference altitudes and 

the wafer top are calculated as above whenever they fall below the cloud 

bottom. However, horizontal coordinates are calculated as a function of 

the wind field only for the wafer top and bottom and any intermediate 

wafer radii are assumed to be defined with respect to the straight line } 
which joins the center of the wafer bottom and a point on the cloud axis 
at the elevation of the lowest intermediate reference plane still in the 
cloud. This process continues for all time steps to stabilization and 
results in a file of descriptions for each wafer in each particle size 
class for all time steps. For the actual calculation of dust density, 
the wafer descriptions are computed for a specific time by an interpola- 
tion process which allows the computations of wafer radius and center as 
a function of altitude. For the calculation of mass concentration at a 
point within a wafer, it is assumed that the mass distribution is uniform 
in the horizontal direction, but exists vertically in equal amounts in 
equal altitude increments. For calculating dust density at a time 

after stabilization, the description of each wafer as it exists at 
stabilization time is used to convert the individual wafer into a stack 
of smaller cylindrical wafers, which are transported individually to the 
time of interest. For the post-stabilization case, the cvlindrical 
wafers are assumed to have their mass uniformly distributed and are 
transported without distortion, 


The user is allowed two methods for specifying the coordinates in 
space at which density calculations will be performed. The first method 
is to simulate a path through space by supplying, as input data to the 
code, a list of coordinates with accompanying times of interest. The 
output for this option is a printed list of the mass densities calculated 
for the separate points at the applicable times. Figure 2 demonstrates 
the computer-generated output for a low yield surface detonation giving 
the mass concentrations at the listed coordinates in space with the 
accompanying times all the same. The mass concentrations are in 10 8 g 
em 3 units, coordinates are in metres, and time in seconds after detona- 
tion. 


The second method provides a cross-sectional view of a cloud at any 
specified time. The user specifies the orientation of a plane through 
ground zero, perpendicular to the ground, and the outer boundaries of a 
rectangle of interest in the plane. Figure 3 illustrates this output 
Option. The detonation was of low yield and a single direction wind 
field was applied in which the velocity increased gradually with altitude. 
The time is 350 seconds, which is before stabilization. Contours of 
equal mass concentrations were drawn directly on the map. Mass con- 
centrations are in units of 10-® g cm 3,altitude and downwind range in 
metres. 
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Figure 2. Dust Concentration at Target Coord 
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IV. USER'S GUIDE TO CASSANDRA 
A. Control Cards for UNIVAC 1108 EXEC 8 


The following sequence of control cards are input for a CASSANDRA 
run on UNIVAC 1108 EXEC 8: 


@ RUN RUN-ID,ACCT-ID,PROJ-1ID,10,200 
@ ASG,A CASSANDRA. 
@ ASG,T 8.,F//POS/10 
@ XQT CASSANDRA.M3DUST 
Insert data deck here for initial conditions, cloud rise, and 
transport. 
S. @ XQF CASSANDRA.MPDQ 
Insert data deck here for CASSANDRA output. 
6. Blank card following last data card. 
7. @ FIN 


B. Data Deck for Execution of CASSANDRA. M3DUST 


1. Deck 1 - Initial Conditions 


Card Sequence Mnemonic Format Read In 
Card 1, File NUMTAP(I) ,T=1,15 1514 MTHRE 
Assignments 
where: 


NUMTAP (1) =ISIN 
NUMT AP (2) =ISOUT 
NUMTAP (3) =IRISE 
and ISIN=Card Reader, set to 5 
ISOUT=Printer, set to 6 
IRISE=Drum assignment, set to 8, the same 
as the control card assignment. 


Card 2 LGO , LTHRU 212 MTHRE 


where: LGO 


= Start point, punch 1 in column 2. 
LTHRU = 


Stop point, punch 2 in column 4. 
Card 3 DETID(J) ,J=1,12 12A6 LINK 1 
where: DETID(J) = Arbitrary 72 character identification 


table for the initial conditions. It 
will be printed as it appears on the 
card, 
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Card 4 N 15 LINK 1 
where: N = The number of data sets for which initial 
condition calculations are desired. If zero 
or negative, 1 is assumed. 


Card 5S IDISTR 15 LINK 1 


where: IDISTR = Preshot soil size - frequency option 
parameter with values as follows: 


IDISTR PARTICLE SIZE DISTRIBUTION 


Ss] Log Norma] 
73 Power Law 
3 Tabular 
Card 6 NDSTR IS LINK 1 
where: NDSTR = Number of entries in the soil particle 


size-frequency distribution. This has 
a maximum value of 200 and is set to 100 
if not specified. 

The following card order is controlled by the value of IDISTR 


For IDISTR $1 


Card 7 W, HEIGHT , USOIL, DMEAN,. 7F10.3 LINK 1 
SD, DNS ,HEW 
For IDISTR = 2 
Card 7 W,HEIGHT ,USOIL,EXPO, 7F10.3 LINK 1 
CAY , DNS , HEW 
Card 8 thru A(I) ,Y(1) , 1=1,NSEG 2E10.4 DSTBN 
CAY + 8 


For IDISTR = 3 | 


Card 7 W,HEIGHT, USOIL, DNS, 7F10.3 LINK 1 
HEW 
Card 8 thru FMASS(I) ,DIAM(I), 2E12.5 LINK 1 
NDSTR + 8 I=1,NDSTR 
Card NDSTR + 9 DIAM(NDSTR + 1) £12.5 LINK 1 
where: W = Total weapon yield in kt. 
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HEIGHT = Height of burst in metres. May be negative for depth of 
burst. 


USOIL = Indicator for soil type. 


USOIL = 0. for siliceous. 
USOIL = 1. for calcareous. 


DMEAN = Mean diameter of the smallest particle size class to be 
used. 
EXPO = Exponent of the frequency function for the power law 
particle size frequency distribution. 
CAY = Coefficient of the frequency function for the power law 
particle size frequency distribution. 
SD = The standard deviation of the log-normal preshot particle 
size frequency distribution. 
DNS = Fallout particle density (g cm 3). Defaults to 2.6. 
HEW = High explosive yield if non-nuclear, 
A = Diameter (microns) of particles from the power law 
particle size distribution curve. 
Y = The percent finer by weight a particle size from the 
power law particle size distribution curve. 
NSEG = IFIX (CAY). 
FMASS = Fraction of the total particulate mass in the ith particle 
Size class. 
DIAM = Diameter of the largest particle in the ith particle size 


class. DIAM (NDSTR + 1) is the diameter of the smallest 
particle in the last (smallest) particle size class. 


2. Deck 2 + GZ Wind Data 


Card Sequence Mnemonic Format _ Read In 
Card 1 NHODO IS SHWIND 
where: 


NHODO = Number of wind levels in wind data deck. If 
NHODO = O no shot-time winds have been 
specified. 
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If NHODO > 0 


Card 2 ZV(J),VX(J) ,VY(J), 3E13.6 SHWIND 
J=1,NHODO 
where: 
ZV = Height of the jth wind stratum. 
VX = X component of the wind velocity of the jth 
stratum (m/s). 
VY = Y component of the wind velocity of the jth 


stratum (m/s), 


3. Deck 3 = LINK 2 Input 


Card Sequence Mnemonic Format Read In 
Card 1 DNID 12A6 LINK 2 
Card 2 KDI, IRAD,KCLD, 614 LINK 2 
KRX, [PAM, KATM 
Card 3 ZBRSTZ E125 LINK 2 
Card 4 SLDTMP E125 LINK 2 
Card 5 FW EVZeS LINK 2 
Card 6 PHI EL2.5 LINK 2 
Card 7 ATID 12A6 LINK 2 
where: 
DNID = Arbritrary cloud-rise module identification. 
KDI = The number of deposit increments per particle 
size class. 
IRAD = Number of cloud wafer radius subdivisions. 


KCLD = Control Index for CRM debug printout. 


KCLD = 1, print data from CRM. 
KCLD = 0, no debug printout. 
KRX = Control Index for RSXP debug printout 
KRX = 0, no printout, KRX = 1, print debug 
data from RSXP. 
IPAM = Particle activity calculation control 
(always zero), 
KATM = Atmosphere printout switch. 
KATM = 0, no printout from ATMR, 
KATM = 1, print complete atmosphere, 
ZBRSTZ = Altitude of GZ (metres above MSL). 
SLDTMP = Particle solidification temperature (°K), 
FW = Fission yield (kt). 
PHI Fraction of FW used to heat air internally. 
ATID = An arbitrary 72 character identification 
for the atmosphere. 


Card 8 FMT 12A6 ATMR 


where: 


FMT = Atmosphere data are read 
by the format specified on the FMT card. The 
eight values are then scaled by the data on 
the SCALE card and are taken to be in the 
order specified on the Nl1---N8 card. 


Cards 9 & 10 SCALE(I) ,I=1,10+ 7F10.5/3F10.5 ATMR 

Card 11 Nl - N8 814 ATMR 
where: 

Nl -- N8 = Atmosphere data sequence indices. 

Card 12 NPVA 14 ATMR 
where: 


NPVA = Number of atmosphere table entries. 


Card 13 to NPVA+11 AP(I),I=1,8 FMT ATMR 


where: 
Ap = Atmosphere data. The data required are: 


ALT = Altitude of the ith level above MSL(metres). 
ATP = Temperature of the ith atmosphere level(°K). 
PRS = Pressure of the ith atmosphere level(mb),. 

RHZ = Atmospheric air density of the ith atmospheric 


level(kg/m3). 
RLH = Relative humidity of the ith atmospheric 
level (percent), 
ETA = Atmospheric dynamic viscosity of the ith 
atmospheric level (kg/(m-s). Does not have 
to be input. 
GRV = Acceleration due to gravity at the ith 
atmospheric level (cm s *). Need not be input. 
SLM = Mean-free-path of air molecules at the ith 
atmospheric level (m). Need not be input. 


then: 


If the data are in the above order, the sequence card would be 
punched 1, 2, 3, 4, 5, 6, 7, 8; however, if not, the proper number 
sequence would be punched and AP(I) is converted as required. For this 
program PRS or RHZ (but not both) need be input. 


4. Deck 4 = Input for Cloud Rise and Point Density 


Gard Sequence Mnemonic Format Read In 
Card 1 DETID(J), J=1, 12 12A6 RSXP 
where: 


DETID = Arbitrary 72 character identification for 
the cloud rise calculations. 


Card 2 IC(J) ,J=1,18 1814 RSXP 


where: 
IC = Control indices. In the current version of 
CASSANDRA only two of the 18 elements of this 
control array are in use. These are IC(5) 
and IC(6). If IC(5)# 0, skips around writing 
wafer history arrays on IRISE tape. If IC(6) # 
0, debug prints complete wafer history. 


Card 3 XGZ,YGZ,TGZ 3E12.5 RSXP 
where: 
XGZ = X - coordinate of ground zero (m). 
YGZ = Y - coordinate of ground zero (m) 
TGZ = Time of detonation (s), 
Card 4 NSUBD,NONE,NTWO,LWONE, 514 RSXP 
LWTWO : 
where: 
NSUBD = Below cloud wafer subdivisions. Default 
value = S. 
NONE = Starting particle size increment. Default 
value = 1. 


NTWO = Ending particle size increment. Default 
value = NDSTR. 

LWONE= Starting wafer bounds. Default value = 1. 

LWTWO= Ending wafer bounds. Default value = KDPST 
where KDPST is set depending upon KDI. If 
KDI > O, KDPST = KDI. If KDI = 0, 
KDPST = 10. 


This is the end of the data decks required for execution of CASSANDRA. 
M3DUST. 


C. Data Deck for Execution of CASSANDRA.MPDQ 


The following input data deck is placed following the @XQT CASSANDRA. 
MPDQ control card. 
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Card Sequence Mnemonic Format Read In 
Card 1 IPDQ 12 PDQXQT 
where: 


IPDQ = Logical unit of Cloud Tape (IRISE). 
Care 2 OPTION 15 PDQXQT 


where: ve 
OPTION = Specifies target option to determine type 
| of calculation and times to be used. 


OPTION = O program stops. 

OPTION = 1 gives tabular output of dust 
loading at specified target 
coordinates. 

OPTION = 2 gives map output of mass 
concentration at points in space. 

OPTION = -1 or #2 are used to repeat 
calculation for a new set of 
coordinates and time. 


IF OPTION = 1 the following cards are input: 


Cards 3 and 4 NTARP,TX(J) ,TY(J) ; I5/4E15.5 PDQXQT 
TZ(J),TRAT(J) ,J=1, 
NTARP 
where: 
NTARP = Number of target coordinates for which 


calculations are to be made. 


TX = X - coordinate relative to GZ (m). 
TY = Y - coordinate relative to GZ (m). 
TZ = Altitude above ground level (m). 


TRAT =Time after detonation at which calculation 
is made (s). 


Card 5 COMBIN, ISTART, ISTOP Slo PDQXQT 


where: 


COMBIN = Control number of particle size classes 
being summed at each target point. If 0, 
indicates summation of all particle size 


classes. 


ISTART = Control number for beginning particle 
size class to be treated. If 0, 
ISTART = NONE. 

ISTOP = Control number for ending particle size 


class to be treated. If 0, ISTOP = NDSTR, 


Card 6 If blank, program terminates since OPTION = 0. 
If -1, program repeats if a new data deck for OPTION 1 follows 
this card. 
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IF OPTION = 2, the following cards are input: 


Card 3 NSTRIP ,NROW, NANGL 315 PDQXQT 
where: 

NSTRIP = Control integer for number of map strips 
to be printed. 

NROW = Control integer for number of rows of 
point density data to be printed. 

NANGL = Control integer for angle relative to 
east of GZ. 


For example, if NSTRIP = 1, NROW = 2 and NANGL = O, a map will be 
produced one page wide, with 15 points in Xdirection and 30 points in 
the Z direction and the angle will be O° from east. 


Card 4 RMIN, RMAX, ALTMIN 6E12.4 PDQXQT 
ALTMAX, DELR, DELALT 


where: 
RMIN = Minimum distance relative GZ (m). 
RMAX = Maximum distance relative GZ (m). 
ALTMIN = Minimum altitude above ground (m), 
ALTMAX = Maximum altitude above ground (m), 
DELR = Horizontal distance between target points (m). 
DELALT = Vertical distance between target points (m). 
One of the above other than RMIN and ALTMIN must 


be input. 
Card 5 TRAT E12.4 PDQXQT 
where: 
TRAT = Time (s) for matrix target representation, 
seconds after burst. 
Card 6 COMBIN, ISTART, ISTOP 815 PDQXQT 
= Same as previously discussed 
Card 7 OPTION I5 PDQXQT 
where: OPTION = 0 program stops. 


OPTION = -2 program repeats for a new time update 


request (Card 5). 
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V. EVALUATION OF CASSANDRA 
A. Dial Pack Simulation 


The validity of the CASSANDRA code was tested by running a 
simulation of the dust cloud resulting from the high explosive event 
Dial Pack. The Dial Pack event was chosen because there was a project 
for the specific purpose of sampling the dust cloud by aircraft traversal 
at a number of times and altitudes. The sampling project resulted in 
calculated values for average mass concentration at various altitudes 
within the cloud. 


In order to simulate the cloud rise dynamics of a high explosive 
(HE) test with a program specifically designed for modelling nuclear 
clouds, it was necessary to modify the assumed energy partition. For 
HE, all of the energy was assumed to heat the cloud. Figure 4 shows 
the cloud top altitude as a function of time, calculated by CASSANDRA, 
together with some observed values for the Dial Pack cloud. Considering 
the irregular nature of the Dial Pack cloud, this is considered to be 
in fair agreement. 


The particle size distribution used to characterize the preshot 
S0il is shown in Figure 5. It is a continuous piecewise power law dis- 
tribution and was constructed to fit the data points shown. The crosses 
represent data obtained from the analysis of a core sample of preshot 
soil. This particular core sample was chosen because the particle size 
range compared favorably with that found in the samples from the air- 
craft traversals and from ground based collectors. The size range used 
was from 0.012 to 2000 microns. 


Due to the lack of shot time wind information for Dial Pack and 
the fact that the average wind changed direction by about 180° shortly 
before detonation, the simulation was run with zero wind. Since the 
figures for comparison are averages over a cloud traversal at a fixed 
altitude, this is not considered a serious drawback. 


2. Hyman, David S., et al, "Dust Cloud Analysis for Event Dial Pack", 


DASA 2694, Headquarters, Defense Atomic Support Agency, Washington, 
DC 20305, May 71. 
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Figure 4, Cloud Rise 
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Figure 5, Dial Pack Particle Size Distribution 
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B. Sample Data Deck for Dial Pack Simulation 


QASGeA CASSANDRA+ 
AaASGeT BeF//POS/10 


axOT CASSANDRAeMSDUST 
5 6 aR 15 16 20 21 22 17 18 
12 
2 
100 
1.9 4.1 200 
2000.0 109.9 
215-0 95.9 
6929 9.5 
1.9 1.6 
in 
459.892 
1473.9 
0.5 
1.0 
NTAL PACK ATMASPHERTC DATA 
(5F10-0/3°2.0) 
1 2 3 4 5 6 7 8 
15 
299,2 931.4 
7A4.555 286.9 650. 
1552.344 281.1 775. 
2326.54 277.5 756 
312.47 27161 638. 
391AR.20 26767 577. 
4722.57 262.5 520. 
5555.59 256.6 4666 
6380.38 250.5 417. 
7246.92 249.5 370. 
AngI5.18 236.5 333. 
9045.95 230.8 258. 
ga4a,31 230.0 253. 
19762018 2ulel 2206 
11712.55 221.1 190. 
ay¥OT CASS*NDPA,MPNOD 
Q 
1 
S 
100 
1500 
2000 
2500 
3000 
10 
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11 12 NUMTAP 
LeO0 


DETID 
N 


IntTsTe? 
NOSTR 
WHUE 

ACL)YO1) 

A(2)Y(2) 

A(3)Y(3) 

A(4)Y(4) 


H4en 


FT! 
SCALF1 
SCALE? 

SEN 

NPVA 

Aqg™ 

AT” 

ATgT™ 

ag 

aT 

ATgT 

art 

AT 

Ag’ 

ATVI 

ATM11 

AT”12 

AT™13 

AT14 

ATYIS 

NETID 

Ic 

X67 evG72"0TSG 

NSB 


Bre 
526 
44. 
Bhe 
39. 
3n. 
39. 
306 
3. 
306 
3ne 
30, 
3M 
30. 
30. 


ODANON SEWN 


IPne 
APTTON 
NTARP 
e 120. 
e 1200. 
° 1200. 
. 120. 
. 1200. 
COMRIN 


CG. Sample Computer Output for Dial Pack Simulation 
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D. Dial Pack Simulation Results 


For Dial Pack, the dust samples were collected by simply opening 
a collector and flying through the cloud. The sample was then weighed 
and a value for average mass concentration obtained by calculating 
the volume swept out by the samples while it was open and dividing the 
mass of the sample by the volume. The dust concentration values from 
sampling which were chosen for comparison are those obtained with an 
aircraft path through the center of the cloud as closely as possible. 
Many of the samples were near the cloud edge and the values fluctuated 
considerably. 


Figure 6 shows the CASSANDRA calculation for the Dial Pack cloud 
at 615 seconds. The aircraft traversal which took place at that time 
was at 2694 metres altitude. The average mass concentration calculated 
from the aircraft sample is 1.6 x 1078 g cm~3, shown in parentheses 
next to the arrow representing the aircraft sampling altitude. The 
value above the aircraft mass concentration is the average CASSANDRA 
calculated value at that altitude. The 2.8 below is the CASSANDRA 
calculated average concentration at 2590 metres altitude. The air- 
craft path was somewhat more toward the cloud edge than the rest so 
the value is lower than the CASSANDRA value. Figure 7 shows the 
CASSANDRA calculation for Dial Pack at 804 seconds. The aircraft 
traversed the cloud at that time at 2542 metres altitude. The average 
mass concentration represented by the aircraft sample was 2.0 x 1078, cm” 
as indicated in parentheses at the arrow representing the aircraft 
position, The average mass concentrations calculated by CASSANDRA 
were 2.0 x 10-8 g cm™? for the altitude of 2590 metres and 
3.0 x 10-8 g cm-? for 2485 metres. These are shown above and below 
the aircraft value. 


3 


CASSANDRA calculations for Dial Pack cloud at 950 seconds, 1304 seconds, 
and 1500 seconds are presented as Figures 8, 9, and 10 respective'y. 
Aircraft dust samples were obtained at these times at altitudes of 
2390 metres, 2237 metres, and 2085 metres respectively. As before, the 
average mass concentrations obtained from aircraft samples are shown 
within parentheses next to the arrow in each figure and CASSANDRA 
average concentrations are printed for altitudes above and below 
aircraft results. 


In conclusion, the CASSANDRA code has been evaluated by comparison 
with available field test data. Tne Dial Pack simulation using 
CASSANDRA has shown excellent agreement with dust concentrations obtained 
by aircraft traversing near or at the center of the Dial Pack dust cloud. 
The agreement constitutes a validation of CASSANDRA for calculating 
dust concentrations. 
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Dial Pack at 950 Seconds Dust Concentration 


Figure 8, 
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Figure 9, 
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Figure 10. 
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